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The two-neutrino double beta decay matrix elements and half- lives of 48 Ca, are calculated within 
a shell-model approach for transitions to the ground state and to the 2 + first excited state of 48 Ti. 
We use the full pf model space and the GXPF1A interaction, which was recently proposed to 
describe the spectroscopic properties of the nuclei in the nuclear mass region A=47-66. Our results 
are T 1/2 (0 + -> 0+) = 3.3 x 10 19 yr and T 1/2 (0 + -> 2+) = 8.5 x 10 23 yr. The result for the decay to 
the 48 Ti + ground state is in good agreement with experiment. The half-life for the decay to the 
2 + state is two orders of magnitude larger than obtained previously. 

£3 ' PACS numbers: 23.40.Bw, 21.60.Cs, 23.40.Hc 
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At present, the double-beta (/?/?) decay is the most sensitive process for direct measurements of the electron neutrino 



presen 
£\j mass. 

13- [1 For der iving the neutrino mass one needs, on one hand, experimental half-lives for the neutrinoless (3(3 
(Qv(3(3) decay mode and, on the other hand, theoretical values of the nuclear matrix elements (NME) entering these 
half-lives formulae. 

After many years of intense investigations and debate on different nuclear structure methods, accurate calculation 
of the NME relevant for (3/3 decay remains a challenging issue. Since many /3j3 emitters are nuclei with open shells, 



the proton-neutron randornphase approximation (pnQRPA) and its extensions, have been the most used methods to 
perform such calculations, [|f- [HI However, due to the significant progress in shell- model (SM) configuration mixing 



(SM 

approaches, there are now calculations performed with these methods for several nuclei. [ 1 lj] - [15l] In spite of their success 
t*"*. . in getting agreement with the experimental half-lives of the two-neutrino (3/3 {2v(3(3) decay mode, both pnQRPA- and 
SM-based approaches have some shortcomings that limit their predictive power for the NME in the case of the 
more interesting 0v(3(3 mode. For example, within pnQRPA methods the NME exhibit a high sensitivity to the 
■ renormalization of the particle-particle strength in the I + channel, while within SM one has to severely truncate the 
model space in order to make the diagonalization procedure tractable. In order to better understand the source of 
uncertainties of the NME calculations for (3/3 decay, a systematic comparison between calculations performed with 
both types of methods is needed. This comparison will become more feasible as the computational power of the shell 
model methods expand to treat larger model spaces. The effective two-body interaction employed is also important, 
since the B(GT + ) strengths are especially sensitive to these interactions. 

SM calculations for (3/3 decay can now be carried out rather accurately for 48 Ca. Zhao, Brown and Richter [TlT ] 
calculated the 2v(3(3 NME of 48 Ca in a large basis SM space using the OXBASH code with the MH (Muto and Horie) 
[r3 | and MSOBEP [l3| two-body interactions. Their predicted Ty 2 is smaller than the experimental one. They also 

made an analysis of the distribution of the B(GT~), B(GT + ) and M^V T components over the Xt excitation energies in 
the intermediate nucleus ( 48 Sc), which helps better understand the quenching of the NME for the 2v(3(3 decay mode. 
Caurier, Poves and Zuker [l4| performed a full pf shell calculation of the NME for the 2v(3(3 decay mode, both for 
the transitions to the g.s. and to the 2+ of 48 Ti. Their calculations were carried out with the ANTOINE code. 19] 
As an effective interaction they used the Kuo-Brown G-matrix \YI\ with minimal monopole modifications, KB3. 18j 
We will discuss their results together with our new results below. 

In this paper we use the recently proposed GXPF1A two-body effective interaction, which has been successfully 
tested for the pf shell ^3|-[23|, to perform 2v/3/3 decay calculations for 48 Ca. Our goal is to obtain the values of 
the NME for this decay mode, both for transitions to the g.s. and to the 2f state of 48 Ti, with increased degree 
of confidence, which will allow us in the next future to address similar calculations for the 0i/(3(3 decay mode of this 
nucleus. [23| The 2v(3/3 transitions to excited states have longer half- lives, as compared with the transitions to the 
g.s., due to the reduced values of the corresponding phase spaces. Positive results for the 2v(3(3 decay of 100 Mo, were 
recently reported. [l9j . 

For the 2v(3(3 decay mode the relevant NME are of Gamow- Teller type, and have the following expressions [lj]- 01 : 
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FIG. 1: The 2v(3(3 decay running matrix element of Eq. (1) as a function of the excitation energy of the intermediate 1 + states 
in 48 Sc, for the GXPF1 and GXPF1A interactions. 
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for the g.s. to g.s. transition, and 

M GT (2 ) - ^ _____ (2) 

for the g.s. to 2^ transition. Here Ek is the excitation energy of the state of 48 Sc and Eq = \Qpfs(Q + ) + AM, 
E-2 = 5<3/3/3(2 + ) + AM. Q/3{3(0 + ) and Q/3/3(2 + ) are the Q-values corresponding to the (3(3 decays to the g.s. and the 
2\ excited state of the parent nucleus ( 48 Ti) and AM is the 48 Ca - 48 Sc mass difference. 
The (3(3 half-live expression is given by 



= Ff \Ml? T (J)f (3) 

where Fj u are the phase space factors[l]: 1.044 x 10~ 17 yr~ x MeV 2 , corresponding to the g.s. to g.s. transition (J=0) 
and 1.958 x 1CP 19 yr~ x MeV 6 , corresponding to g.s. to 2\ transitions (J=2), respectively. 

The calculations were carried out in the full pf model space using the CMISHSM shell model-code [2^] and the 
GXPF1A interaction. The most recent effective Hamiltonians, GXPF1 [^-[IH and GXPF1A H2 are derived from 
a microscopic calculation by Hjorth- Jensen based on renormalized G matrix theory with the Bonn-C interaction [25| . 
and are refined by a systematic fitting of the important linear combinations of two-body matrix elements to low-lying 
states in nuclei from A=47 to A=66. GXFP1A addresses some shortcoming of the GXPF1 interaction for the region 
of the neutron rich Sc, Ti and Ca isotopes[22] that is relevant for this study. An advantage of using the full pf model 
space is that the Ikeda sum rule is exactly satisfied. 

In the calculation of the NME, Eqs. (1-2), we used the standard quenching factor of 0.77 for the <tt operator. [3| 
We used up to 250 intermediate 1 + 48 Sc states in the sum. They exhaust nearly the entire B(GT) sum rules for the 
transitions from 48 Ti and 48 Ca: 1.59 out of the exact 1.6 for Ti and 22 out the exact 24 for Ca. 

We also tested the validity of the quenching factor of 0.77 by comparing the beta decay probabilities for the 48 
Sc(6 + ) — s- 48 Ti(6 + ) transitions with the experimental data. [3(| The results are presented in Table I and confirm that 
this value reasonably describes the B(GT) quenching for this mass region. 

The running NME, Mq'j, of Eqs. (1-2), as a function of the excitation energy of the 1 + states in 48 Sc are presented 
in Figs. 1 and 2, respectively. The convergence trends are similar to the ones found in Refs. [ll|, H3' and it is also 
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FIG. 2: The 2vf3f3 decay running matrix element of Eq. (2) (up to a factor of \/5) as a function of the excitation energy of the 
intermediate 1 + states in 48 Sc, for the GXPF1 and GXPF1A interactions. 
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FIG. 3: The 2u/3/3 decay running matrix element of Eq. (1) as a function of the excitation energy of the intermediate 1 + states 
in 48 Sc compared with the similar sum where the absolute values of the GT matrix elements are used. 



supported by the nearly exhausted sum rules. It is also clear that the phases of the intermediate states in the double 
sum play an essential role: the contribution of the intermediate states is not coherent and the sum is not continuously 
increasing. [5^] This is illustrated more clearly in Fig. 3 where the coherent sum (signed amplitudes) is compared 
to the incoherent sum (absolute GT amplitudes). Fig. 4 presents the relevant B(GT) transition probabilities. We 
note that most of the positive contribution to the double-beta matrix element to the 48 Ti ground state shown in 
Fig. 1 comes from two out of the five intermediate 1 + states below 5 MeV excitation in 48 Sc. For the product 



TABLE I: Theoretical end experimental log ft for the Sc^ Ti beta decay transitions. 
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3.333 
3.508 


5.247 
6.083 


5.532 
6.010 
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[(0/||CTT-||l£)][(l£||£7T-||0i)] the lowest 1+ (at 2.5 MeV) gives [0.185] [1.15] = 0.122 and the third 1+ (at 3.8 MeV) 
gives [0.42] [0.35] = 0.147. 

Recent experiments have attempted to extract B(GT) values from 48 Ca( 3 He,t) 48 Sc and 48 Ti(d, 2 He) 48 Sc reaction 
cross sections (see e.g. Fig. 4 of a recent review, Ref. |3l|). The results obtained for the lowest strong 1 + observed in 
48 Ca( 3 He,t) at 2.5 MeV [B(GT) 1 / 2 for the above product] is |0.12(3)| x |0.95(5)| = |0.11(3)|, in good_agreement with 
theory given the uncertainties that exist in extracting B(GT) from charge-exchange cross sections [32] (the state at 
2.2 MeV in 48 Ti(d, 2 He) associated with 1 + does not have a correspondence in the theory - it is near a state previously 
assigned 3 + in the literature and its J 77 value should be confirmed.) The double-beta strength associated with the 
theoretical state at 3.8 MeV appears to be spread over several states near 3 MeV in experiment. 

Using the results from Figs. 1-2 one gets the following converged results for the the 2v(3f3 matrix elements: 

• | M(0+ -> 0+) | = 0.0539 MeV" 1 for GXPF1A, and | M(0+ -> 0+) | = 0.0635 MeV" 1 for GXPF1; 

• | M(0+ -> 2+) | = 0.0122 MeV" 3 for GXPF1A, and | M(0+ -> 2+) | = 0.0129 MeV" 3 for GXPF1. 

Using these matrix elements and the phase factors of Ref. [l| in Eqs. (3) one gets for the 2v(3j3 decay half-lives: 

• T 1/2 (0+ -> 0+) = 3.3 x 10 19 yr for GXPF1A, and T 1/2 (0+ 0+) = 2.4 x 10 19 yr for GXPF1; 

• T 1/2 (0+ 2+) = 8.5 x 10 23 yr for GXPF1A, and T 1/2 (0+ 2+) = 7.5 x 10 23 yr GXPF1. 

Our value for the T 1/2 (0+ -> 0+), corresponding to the | M(0+ -»• 0+) | = 0.0539 MeV" 1 NME calculated with 
GXPF1A interaction, is within the experimental range [29]: (4.3^3'g) x 10 19 yr. The calculations performed with the 
GXPF1 interaction seem to give a larger value for the NME that leads to a half-life value which isiust below the 
present experimental range. Comparing our results to the previous similar ones of Refs. [HI and [l| we note that 
the calculations of Zhao, Brown and Richter are performed in a restricted pf model space, and they found a NME 
of 0.07 MeV" 1 as their best value. Furthermore, using a phase space factor slightly different from ours, they obtain 
^i/2(0 + ~~ * + ) = 1.9 x 10 19 yr. This half-life is about half of our value and is significantly below the experimental 
range. 

Caurier, Poves and Zuker found | M(0 + — > + ) | = 0.0402 MeV -1 in their work [lj]. We repeated their calculations 
as described in Ref. [l4j, and obtained 0.047 MeV" 1 for the same NME, value that is in agreement with the NME 
reported by Nowacki in Ref. [l5j |. and which differs by about 13% from ours. This value is also within the present 
experimental range. 

For the g.s. to the 2+ transition we obtained a NME value, which is about half the numerical value reported by 
Caurier, Poves and Zuker in Ref. Q. However, Ref. [3 used the same Eq. (1) for the | M[0+ -> 2+) | NME, 
instead of our Eq. (2), which is recommended in the literature 0, Based on Eq. (1), Ref. [lj] suggests that the 
+ — > 2 + decay rate is about 3% of the + — > + decay rate. Our value obtained for the 2v(3(3 half-life corresponding 
to the transition to the 2\ excited state of 48 Ti from Eq. (2) is about four orders of magnitude larger than that for 
the g.s. to g.s transition. 

In conclusion we calculated the NME and half-lives for 2vf3j3 decay of 48 Ca within a SM approach in the full pf 
model space. We calculated both the g.s. to g.s. and g.s. to 2]*~ excited state transitions. We use for the first time[33j] 
in such calculations the two versions of GXPF1 two-body interaction, which were recently proposed and successfully 
used to reproduce the spectroscopic properties of many nuclei in the nuclear mass range A=47-66. Our results are 
based on 250 1 + intermediate states in 48 Sc nucleus which are enough to exhaust almost the entire B(GT) sum rules 
for the transitions from 48 Ti and 48 Ca. We also checked the validity of 0.77 quenching factor for the Gamow- Teller 
operator used in the SM calculations by comparing the calculated beta transitions 48 Sc^ 48 Ti with the experimental 
ones. The best values, we propose for NME are | M(0+ -> 0+) | = 0.0539 MeV" 1 and | M(0+ 2+) | = 0.0122 
MeV -3 , which were obtained using the GXPF1A interaction. They correspond to Ti/ 2 (0 + — > + ) = 3.3 x 10 19 yr 
and T 1/2 (0+ 2+) = 8.5 x 10 23 yr, respectively. Future experiments on (3/3 decay of 48 Ca, CANDLES [H[ and 
CARVEL (ifl, may reach the required sensitivity of measuring such transitions and our results could be useful for the 
planning of these experiments. 

MH and BAB acknowledge support from the NSF grant PHY-0555366. SS acknowledge travel support from the NSF 
grant INT-0070789. 
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FIG. 4: The relevant B(GT) probabilities for transitions to the 1 + states in 48 Sc. 
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